Asymptotic step profiles from a nonlinear growth equation for vicinal surfaces 
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We study a recently proposed nonlinear evolution equation describing the collective step meander 
on a vicinal surface subject to the Bales-Zangwill growth instability [0. Pierre- Louis et at, Phys. 
Rev. Lett. 80, 4221 (1998)]. A careful numerical analysis shows that the dynamically selected step 
profile consists of sloped segments, given by an inverse error function and steepening as \ft, which are 
matched to pieces of a stationary (time-independent) solution describing the maxima and minima. 
The effect of smoothening by step edge diffusion is included heuristically, and a one-parameter 
family of evolution equations is introduced which contains relaxation by step edge diffusion and by 
attachment-detachment as special cases. The question of the persistence of an initially imposed 
meander wavelength is investigated in relation to recent experiments. 
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I. INTRODUCTION 

Ten years ago, Bales and Zangwill |l|] predicted that a 
growing vicinal surface should undergo a step meander- 
ing instability when kinetic step edge barriers suppress 
the attachment of atoms to descending steps ||. The 
instability has meanwhile been observed in experiments 
HU and Monte Carlo simulations ||, and a number of 
theoretical studies have been devoted to the nonlinear 
evolution of the surface both in the presence || and ab- 
sence HQ of desorption || . 

Since linear stability analysis shows the in-phase mode 
of the collective step meander to be the most unstable 
Q , the two-dimensional surface morphology can be rep- 
resented by a one-dimensional function £(x, t) describing 
the displacement of the common step profile from the 
flat straight reference configuration £ = 0, with the x- 
axis oriented along the step For the case of infinite 
step edge barriers, attachment-detachment kinetics and 
no desorption, the nonlinear evolution equation 
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increase of the meander amplitude as \ft at fixed wave- 
length X u , as well as the formation of spike singularities 
at maxima and minima of The latter is surprising be- 
cause the second term on the right hand side of (0) would 
be expected to suppress such rapid variations of the step 
curvature. 

Here we revisit the problem using a more accurate 
numerical algorithm We demonstrate that the 

step profile remains smooth near maxima and minima, 
where it approaches asymptotically a stationary (time- 
independent) solution of ([!]), while the sides of the profile 
follow a separable solution with an amplitude of order \ft. 
The matching of the two solutions occurs near the point 
of maximum slope. We further show heuristically how the 
effect of step edge diffusion can be included in the the- 
ory, and introduce a generalized evolution equation which 
contains edge diffusion and attachment-detachment ki- 
netics as special cases. Finally, we address the question 
to what extent an initially imposed meander wavelength 
different from X u is preserved under the time evolution. 
This is relevant in view of the recent experiments of 
Maroutian et al. B. 



(1) 



was proposed in Ref. |?J (subscripts denote derivatives). 
It can be derived from the Burton-Cabrera-Frank (BCF) 
theory of growth on vicinal surfaces [ I1J using a singular 
multiscale expansion Jt],[i^] in e 1 / 2 , where e = VtFi 2 / D is 
the Peclet number. Here F is the deposition flux, D the 
in-plane surface diffusion coefficient, £ the nominal step 
spacing, and the atomic area. The coefficients in Eq. (|l ) 
are given by a = QF£ 2 /2 and (3 = VL 2 D£<yc a Jk B T, with 

7 and c cq referring to the step stiffness and the equilib- for a classical particle of mass (3 moving in the potential 



II. SHAPE SELECTION 

Before presenting the numerical results we recapitulate 
the two classes of analytic solutions to ([!]) which were 
found in Qj. Stationary solutions are obtained by set- 
ting the mass current along the step (the quantity inside 
to curly brackets on the right hand side of (Q)) to zero. 
In terms of m(x) — C, x / yl + Cx the stationarity condi- 
tion reduces to Newton's equation (3d 2 m/dx 2 = —dU/dr 



rium adatom density, respectively. 

According to (|l|), the straight step is linearly unstable 
against perturbations with wavelengths larger than A c = 
2ny / f3/a, with a fastest growing wavelength A M = V2X C . 
To explore the nonlinear regime, in Ref. a numerical 
integration of Eq.(nh was carried out which showed an 
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which can be solved by quadra- 



tures. One thus obtains a one-parameter family of peri- 
odic profiles Cs( x ) which are most conveniently parame- 
terized by the maximum slope S = max x Q x , and which 
have been described previously in the context of a differ- 
ent surface evolution equation The amplitude A(S) 
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FIG. 1. The evolution of the step profile starting from 
a flat initial condition with small random fluctuations. The 
upper figure shows the case of attachment-detachment kinet- 
ics (Eq. (|l])) at times t = 36,64,110,183, the lower figure 
the case of edge diffusion (Eq.(tf) with n — 1/2) at times 
t = 20,60,112,200. Subsequent profiles have been shifted 
in the (^-direction. I n all figures spatial variables have been 
scaled by \ c /2ir = fi/a and time by fi/a 2 . 

is an increasing function of S, while the wavelength A(5 I ) 

decreases with increasing S, starting out at A(0) = A c . 

For S —> oo finite limiting values A(oo) = 

A(oo) = y/2nf3/a r(3/4)/r(5/4) » 0. 5393527.. A c are 

approached. 

The separable solution of interest reads fill"! 

C(x,t) = 2y/rtexr 1 (l-4\x\/\ a ), (2) 

-A s /2 < x < X s /2, where eri(z) = (2/Vtt) / * dy e"^, 
and the wavelength A s is arbitrary. Equation ^ solves 
(|I|) exactly in the limit t — > oo, when the second term 
on the right hand side becomes negligible compared to 
the first, and the evolution equation reduces to £ t — 
— {a/C,x)x- The solution (Q) is singular near the maxima 
and minima, where it diverges as £ ~ ±^/ln(l/|a; — xoQ, 
x = 0,±X s /2. 

In Figure [l] we show results of a numerical solution of 
(Q), starting from a small amplitude random initial con- 
dition. To secure good numerical stability we used a fully 
implicit, backwards Euler algorithm for integration. The 
algorithm was implemented on an adaptive grid in or- 
der to obtain sufficient lateral resolution at the singular 
points. 

A regular meander pattern of wavelength A„ develops, 
with an amplitude growing indefinitely as \ft. Closer 
inspection reveals that the sides of the profile follow the 
separable solution (Figure |2|) , while near the maxima and 
minima smooth caps appear which approach pieces of the 
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FIG. 2. The asymptotic form of the scaled profile C/Vi f° r 
Eq- (Q) (long dashes) and Eq.(Q) with n = 1/2 (short dashes). 
Full line is the separable solution (0) with A 3 equal to the 
total meander wavelength A = 18y f3/a w 2\ u . 



stationary solutions (Figure ||). This can be understood 
by noting that the mass current along the sides of the 
profile vanishes as l/y/i according to (0), and therefore 
the stationarity condition is asymptotically satisfied; we 
have checked that the deviation from the stationary pro- 
file which is discernible in Figure || vanishes as \j\ft. 
Since the slope of (||) increases monotonically upon ap- 
proaching an extremum while it decreases for the sta- 
tionary profiles, the matching of the two solutions occurs 
near the point of maximum slope. For t — > oo the slope 
of the separable solution diverges, hence the cap profile 
approaches the limiting stationary solution Coo(a^), and 
the length of the cap becomes A(oo)/2. The rescaled step 
profile £(&, t)/s/i approaches an invariant shape in which 
the cap appears as a flat facet. The wavelength X s of the 
separable solution depends on the cap length and on the 
total meander wavelength A, and is fixed by mass balance 
requirements Q ; for large total wavelength A s — * A (see 
Figure I). 
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FIG. 3. The form of the caps for Eq.(§) 
at t = 500, 1500, 2500, 3400. The dashed lines are the sta- 
tionary solutions corresponding to the maximum slope in the 
profile. 
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FIG. 4. Asymptotic form of the step profile for n — 1/2 
(upper) and n = (lower). Full line is the profile and dashed 
line the slope £ x ■ Circles represent the stationary solution and 
diamonds the corresponding slope. The separable solution, 
still present in the sloped regions of the profile for n = 1/2, 
has vanished for n = 0. 



III. STEP EDGE DIFFUSION AND A 
GENERALIZED EVOLUTION EQUATION 

On many fee metal surfaces, diffusion along step edges 
is the fastest kinetic process which therefore provides the 
dominant step smoothening mechanism To see how 
Eq.(|l]) has to be modified to take this effect into account, 
we note first that the second, relaxational term on the 
right hand side can be rewritten in a geometrically covari- 
ant form as (crfj, s ) x , where /i = tt^/K is the step chemical 
potential Jig], k = —(l + (l)~ 3 / 2 ( xx is the step curvature, 
s = j dx i/lT^ is the arclength along the step, and 
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is a mobility. The ^-dependence in (||) reflects the as- 
sumed relaxation kinetics [m], in which mass exchange 
between different parts of the step occurs through de- 
tachment followed by diffusion over the terrace, reflec- 
tion at the descending step, and re-attachment (case E 
of Jl8| ) . The factor 1 / y^l + ( 2 has a simple geometric 
interpretation |l7j ]: For a deformed in-phase step train 
the distance to the nearest step, measured along the step 
normal, is i/y/l + C 2 rather than i. 

For relaxation through step edge diffusion the mobility 
is clearly independent of the step distance, and is given 
by [|l8| a = D e £lc e /kBT, where D e and c e denote the edge 
diffusion coefficient and the equilibrium concentration of 
edge atoms, respectively. When edge diffusion dominates 



(i.e. a 3> cr), the appropriate nonlinear growth equation 
should thus be given by (|ll) with the second term replaced 
by {d-fj, s ) x = ([1 + C, 2 x Y^ 2 ayL x ) x . This is confirmed by 
the explicit derivation of Gillet et al. , who also stud- 
ied the crossover between attachment-detachment kinet- 
ics and edge diffusion. Here our primary goal is to gain 
further insight into the shape selection mechanism. This 
has lead us to consider the generalized class of equations 
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which reduces to (Q) for n = 1 and describes relax- 
ation through step edge diffusion when n — 1/2 and 
(3 = £l 2 D c c c ^ jhsT '. Below we discuss the properties of 
([|) for general n, keeping in mind that the cases n = 1/2 
and n = 1 are of immediate physical relevance. 

The separable solution (|) becomes exact in a limit 
where the relaxation term in (|l|) can be neglected, hence 
it remains a valid asymptotic solution also of (Q) for 
n > —1/2; for n < —1/2 the relaxation term can 
never be ignored. The stationary solutions of ([|) can 
be analyzed in terms of the same mechanical analogy 
described above, the particle potential being given by 
U(m) = -a(l - m 2 ) 3 / 2 -'7(3 - 2n). For 1/2 < n < 3/2 
the behavior is analogous to that for n = 1: The wave- 
length A(S) is a decreasing function of the maximal slope 
5, and wavelength and amplitude reach finite values 
A(oo) = \/8/3/qV3 - 2n/(2n - 1) and 



A(oo) = V2tt(3 - 2n){j3/a) 



r[(2n+ l)/4] 
r[(2n + 3)/4] 



(5) 



for S — > 00. Thus the asymptotic step profiles look sim- 
ilar to those generated by Eq. ([!]) , with the length of the 
cap decreasing with increasing n. As n — > 3/2 the cap 
length, given by A(oo)/2, vanishes. For n > 3/2 we 
therefore expect true spike singularities to develop at the 
maxima and minima of the profile. Using that the slope 
imposed by the separable solution (||) grows as S ~ Vt, 
we predict that the curvature at the extrema diverges as 

f (2n-3)/4_ 

A numerical solution for the case of edge diffusion 
(n = 1/2) is shown in the lower panel of Figure ffl. For 
n = 1/2 the potential U(m) is harmonic and hence the 
wavelength A(S) = X c independent of S. The ampli- 
tude of the stationary profiles diverges as A(S) ~ lnS*, 
leading to a corresponding increase of the cap height as 
Int. Since this is still small compared to the overall pro- 
file amplitude, the caps nevertheless appear as flat in the 
rescaled shape C/Vt _(Figure ||; a detailed view of the cap 
is shown in Figure ^). This remains true in the entire 
interval —1/2 < n < 1/2, where the cap length re- 
mains finite and the cap height grows as t( 1-2 ™)/ 4 . How- 
ever a qualitative change in the profile evolution occurs 
at the value n c ps 0.2283 where the asymptotic station- 
ary wavelength (||) becomes equal to the most unstable 
wavelength X u , which sets the lateral length scale in the 



3 




0.0 1 ■ ' ■ ' 1 

0.0 10.0 20.0 

X 

FIG. 5. Spontaneous creation of an extra meander period 
for the case of step-edge diffusion (n = 1/2). The initial 
wavelength of the profile is \i = 25W~f3/a > 3A C . 

early stages of growth. For n < n c we expect to see 
an intermediate coarsening regime in which the lateral 
length scale increases from X u to A(oo). Coarsening to 
arbitrarily large length scales, similar to what is observed 
in related evolution equations for one-dimensional unsta- 
ble growth |lS[ ], sets in at n = —1/2, where (^) diverges. 
Throughout the regime n < n c the evolving profile is 
describable in terms of stationary solutions, and the sep- 
arable solution (0) no longer plays any role (Figure ||). 



IV. PERSISTENCE OF THE INITIAL 
WAVELENGTH 

Finally, we address the recent experiments Q on sur- 
faces vicinal to Cu(100), in which the meander wave- 
length was measured as a function of temperature, and 
it was concluded that the observed behavior is inconsis- 
tent with the theoretical prediction for the linearly most 
unstable wavelength A u . Maroutian et al. [Q therefore 
proposed that the meander wavelength is set by the nu- 
cleation length describing the distance between the one- 
dimensional nuclei appearing on a flat step in the early 
stage of growth, which can be considerably larger than 

A necessary consistency requirement for this scenario 
is that an initially imposed meander wavelength \ > X u 
persists during the nonlinear evolution. We have there- 
fore numerically integrated Eqs.([l], |J) starting from a 
sinusoidal initial condition with varying wavelength A^. 
We do find that a range of wavelengths can be preserved 
during growth. This is reasonable in view of the analysis 
presented above, which shows that asymptotic profiles, 
composed of the separable solution (||) and a stationary 
cap, can in principle be constructed for arbitrary wave- 
length (see e.g. Figure ||). However, when A, exceeds A c 
by more than a factor of 3, so that an additional mean- 
der fits between the maxima and minima of the profile, 
the wavelength spontaneously decreases to a value near 
X u (Figure ||) . This result contradicts the assumption of 
Q that initial wavelengths much larger than X u persist, 



but it should not be overemphasized: Clearly processes 
which involve a change in the collective meander wave- 
length may not be accurately described in a model which 
assumes in-phase meandering from the outset. 



V. OUTLOOK 

In conclusion, we have described an unusual shape se- 
lection scenario for a class of physically motivated growth 
equations. A number of issues remain to be clarified. 
Mathematically, the behavior in the region where sepa- 
rable and stationary solutions match needs further inves- 
tigation; our numerical work indicates the appearance of 
singularities in higher derivatives of Also the dynam- 
ics in the singular regime n > 3/2 and in the coarsening 
regime n < n c of Eq. (^) deserves attention. Physically, it 
is imperative that the predictions of the one-dimensional 
equations for the in-phase step meander be confirmed 
by more complete descriptions of the growing surface, as 
provided by two-dimensional continuum equations and 
Monte Carlo models Q , in order to assess their ultimate 
relevance for the experimentally observed morphologies. 
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